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Abstract In this paper, the effect of self-gravity on the protoplanetary discs is investi- 
gated. The mechanisms of angular momentum transport and energy dissipation are as- 
sumed to be the viscosity due to turbulence in the accretion disc. The energy equation 
is considered in situation that the released energy by viscosity dissipation is balanced 
with cooling processes. The viscosity is obtained by equality of dissipation and cool- 
ing functions, and is used for angular momentum equation. The cooling rate of the flow 
is calculated by a prescription, du/dt = —u/tcooI, that u and TcooI are internal energy 
and cooling timescale, respectively. The ratio of local cooling to dynamical timescales 
^Tcooi is assumed as a constant and also as a function of local temperature. The solu- 
tions for protoplanetary discs show that in situation of HtcooI = constant, the disc does 
not show any gravitational instability in small radii for a typically mass accretion rate, 
M ~ 10~^MQyr^^, while by choosing HtcooI as a function of temperature, the gravi- 
tational instability for this amount of mass accretion rate or even less can occur in small 
radii. Also, by study of the viscous parameter a, we find that the strength of turbulence 
in the inner part of self-gravitating protoplanetary discs is very low. These results are 
qualitatively consistent with direct numerical simulations of protoplanetary discs. Also, 
in the case of cooling with temperature dependence, the effect of physical parameters 
on the structure of the disc is investigated. The solutions represent that disc thickness 
and Toomre parameter decrease by adding the ratio of disc mass to central object mass. 
While, the disc thickness and Toomre parameter increase by adding mass accretion rate. 
Furthermore, for typically input parameters such as mass accretion rate lO^^A/0?;r~^, 
the ratio of the specific heats 7 = 5/3, and the ratio of disc mass to central object mass 
q = 0.1, the gravitational instability can occur in whole radii of the discs excluding very 
near to the central object. 
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1 INTRODUCTION 

Accretion discs are important for many astrophysical phenomena, including protoplanetary systems, 
different types of binary stars, binary X-ray sources, quasars, and Active Galactic Nuclei (AGN). 
Historically, theory of accretion discs had concentrated in the case of non self-gravitating and occa- 
sionally the effect of self-gravity had studied (Paczyriski 1978; Kolykhalov & Sunyaev 1979; Lin & 
Pringle 1987, 1990). On the other hand, in recent years, the importance of study of disc self-gravity 
has increased, especially in the protostellar discs and Active Galactic Nuclei (AGN) discs. It can be 
due to increase of computational resources in simulation of self-gravitating accretion discs and the ob- 
servational evidences that have confirmed the existence of self-gravity on all scale discs, from AGN 
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to protostars (Lodato 2008 and references therein). Also, it appears the development of gravitational 
instability is important for cool regions of accreting gas that angular momentum transport by magneto- 
rotational instability (MRI) becomes weak (Fleming et al. 2000; Masada & Sano 2008; Faghei 2011) 
and angular momentum can transport by gravitational instability. 

The structure of self-gravitating discs has been studied both through self-similar solutions assuming 
steady and unsteady state (Mineshige & Umemura 1996, 1997; Tsuribe 1999; Bertin & Lodato 1999, 
2001; Shadmehri & Khajenabi 2006; Abbassi et. al. 2006; Shadmehri 2009) and through direct numer- 
ical simulations (Gammie 2001; Rice et al. 2003, 2005, 2010; Rice & Armitage 2009; Cossins et al. 
2010; Meru&Bate 2011a). 

Mineshige & Umemura (1996) investigated the role of self-gravity on the classical self-similar 
solution of advection dominated accretion flows (ADAF, Narayan & Yi 1994) and found global one- 
dimensional solutions influenced by self-gravity in both the radial and the perpendicular directions of 
the disc. They extended the previous steady state solutions to the time-dependent case while the effect 
of the self-gravity of the disc was taken into account. They used an isothermal equation, and so their 
solutions describe a viscous accretion discs in the slow accretion limit. Tsuribe (1999) studied unsteady 
viscous accretion in self-gravitating discs. Taking into account the growth of the central point mass, 
Tsuribe (1999) derived a series of self-similar solutions for rotating isothermal discs. The solutions 
showed, as a core mass increases, the rotation law changes from flat rotation to Keplerian rotation in 
the inner disc and in addition to the central point mass, the inner disc grows by mass accumulation due 
to the differing mass accretion rates in the inner and outer radii. Bertin & Lodato (1999) considered 
a class of steady-state self-gravitating accretion discs for which efficient cooling mechanisms are as- 
sumed to operate so that the disc is self -regulated at a condition of approximate marginal Jeans stability. 
They investigated the entire parameter space available for such self -regulated accretion discs. In another 
study, Bertin & Lodato (2001) followed the model that, when the disc is sufficiently cold, the stirring 
due to Jeans-related instabilities acts as a source of effective heating. The corresponding reformulation 
of the energy equations, they demonstrated how self-regulation can be established, so that the stability 
parameter Q is maintained close to a threshold value, with weak dependence on radius. Abbassi et al. 
(2006) studied the effect of viscosity on the time evolution of axisymmetric, polytropic self-gravitating 
discs around a new born central object. Thus, they ignored from the gravitational effect of central ob- 
ject and only self-gravity of the disc played an important role. They compared effects of a-viscosity 
prescription (Shakura & Sunyaev) and /3-viscosity prescription (Duschel et al. 2000) on disc structure. 
They found that accretion rate onto the central object for /3-discs more than a-discs at least in the outer 
regions where /3-discs are more efficient. Also, their results showed gravitational instability can occur 
everywhere on the /3-discs and thus they suggested that /3-discs can be a good candidate for the origin 
of planetary systems. Shadmehri & Khejenabi (2006) examined steady self-similar solutions of isother- 
mal self-gravitating discs in the presence of a global magnetic field. Similar to Abbassi et. al. (2006) 
they neglected from the mass of the central object to the disc mass. By study of Toomre parameter they 
showed that magnetic field can be important in gravitational stability of the disc. 

An accretion discs can become gravitationally unstable if Toomre parameter becomes smaller than 
its critical value, Q < Qcrit (Toomre 1964). For axisymmetric instabilities Qcrit ~ 1, while for non- 
axisymmetric instabilities Qcrit values as high as 1.5 — 1.7 (Durisen et al. 2007). One possible outcome 
is that unstable discs fragment to produce bound objects and has been suggested as a possible mech- 
anism for forming giant planets (Boss 1998, 2002). However, recently it has been realized that above 
condition is not sufficient to guarantee fragmentation. Gammie (2001) showed that in addition to the 
above instability criterion, the disc must cool at a fast enough rate. Let the cooling timescale TcooI be 
defined as the gas internal energy divided by the volumetric cooling rate. For power-law equation of 
state and with TcooI prescribed to be some value over a annulus of the disc, the thin shearing box simula- 
tions of Gammie (2001) show that fragmentation occurs if and only if flTcooi ^ l^crit, where Pcrit ^ 3 
and n is angular velocity of the disc or inverse of dynamical timescale Tdyn = The critical value 
of flTcooi can be somewhat larger than three for more massive and physically thicker discs (Rice et 
al. 2003), larger adiabatic index (Rice et al. 2005), and more resolution of simulations (Meru & Bate 
2011b). Cossins et al. (2010) by SPH simulation studied the effects of opacity regimes on the stability 



Self-gravitating protoplanetary discs 



3 



of self-gravitating protoplanetary discs to fragmentation into bound objects. They showed that flvcooi 
has a strong dependence on the local temperature. As, they found that without temperature dependence, 
for radii < 10 AU a very large accretion rate 10~^MQyr^^ is required for fragmentation, but that this 
is reduced to 10^'' with cooling of dependent on temperature. 

As mentioned, typically semi-analytical studies of self-gravitating discs are regarding polytropic 
discs (Abbassi et al. 2006), isothermal discs (Mineshige & Umemura 1996, 1997; Tsuribe 1999; 
Shadmehri & Khajenabi 2006), ADAFs in the extreme limit of no radiative cooling (Shadmehri 2004), 
and discs without central object (Mineshige & Umemura 1996, 1997; Tsuribe 1999; Shadmehri & 
Khajenabi 2006; Abbassi et al. 2006). In this paper, it will be interesting to understand under which con- 
ditions gravitational instability can occur in accretion discs by a suitable energy equation and assuming 
a Newtonian potential of a mass point that stands in the disc centre. Thus, to obtain these conditions, we 
will use a prescription for cooling rate that is introduced by Gammie (2001), du/dt = —u/tcooU that u 
and Tcooi are internal energy and cooling timescale, respectively. The ratio of local cooling to dynamical 
timescales ^TcooI is assumed a power-law function of temperature in adapting Cossins et al. (2010), 
^Tcooi = Po{T/To)^ , where Tq and S are free parameters, and /3o is free parameter in Gammie (2001). 
In 5 = 0, flTcooi reduces to Gammie (2001) model that ^TcooI is a constant, while non-zero S is qual- 
itatively consistent with results of Cossins et al. (2010). We will examine the effects of 6 parameter on 
gravitational stability of disc. We will show that the present model is qualitatively consistent with direct 
numerical simulations (Rice & Armitage 2009; Cossin et al. 2010; Rice et al. 2010) and can provide 
conditions that gravitational instability occur in whole radii excluding very near to the central object. 

In section 2, the basic equations of constructing a model for steady self-gravitating disc will be 
defined. In section 3, we will find asymptotic solutions for outer edge of the disc. In section 4, by exploit 
of asymptotic solutions as boundary conditions for system equations, we will investigate numerically 
the effects of physical parameters on structure and stability of the disc. The summary and discussion of 
the model will appear in section 5. 

2 BASIC EQUATIONS 

We use cyUndrical coordinate (r, ip, z) centered on the accreting object and make the following standard 
assumptions: 

(i) The flow is assumed to be steady and axisymmetric dt = d^p = 0, so all flow variables are a function 
of r and z ; 

(ii) The gravitational force of central object on a fluid element is characterized by the Newtonian poten- 
tial of a point mass, = — GA/* / r, with G representing the gravitational constant and standing 
for the mass of the central star; 

(iii) The equations written in cylindrical coordinates are integrated in the vertical direction, hence all 
quantities of the flow variables will be expressed in terms of cylindrical radius r; 

The governing equations on the self-gravitating accretion disc for such assumptions is as follows. The 
continuity equation is 

Id, 

-— rSv, =0, (1) 
r dr 

where Vr is the radial infall velocity and S is the surface density, which is defined as E = 2ph, and p 
and h are density and the disk half-thickness, respectively. The half-thickness of the disc with assume 
of hydrostatic equilibrium in vertical direction is h = Cg/il, where is sound speed, which is defined 
as ~ p/p, p being the gas pressure and ft represents angular velocity of the flow. The equation (1) 
implies that 

Af = — 27rrSvr = constant 

where M is the mass accretion rate and is a constant in the present model. The simulation results of 
protoplanetary discs show that the disc reaches a quasi-steady state in 20000 years or less and might 
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imply that these systems are rarely out of equilibrium. Also, the simulations show that the mass of the 
disc redistribute itself to produce a state in which the accretion rate. A/, is largely independent of r (Rice 
& Armitage 2009, Rice et al. 2010). Thus, we can use the mass accretion as a constant and it can not be 
a limitation for the present model. The momentum equations are 



dVr 1 d 2n 



G 



A/* + M{r) 



(2) 



T,Vr — {r il) ^ - — 
dr r dr 



vY.r-^ 



dr 



(3) 



where ly is the kinematic viscosity coefficient, and 7 is the adiabatic index, and M{r) is the mass of 
a disc within a radius r. As in Mineshige & Umemura (1997), we adopt the monopole approximation 
for the radial gravitational force due to the self-gravity of the disc, which considerably simplifies the 
calculations and is not expected to introduce any significant error as long as the surface density profile 
is steeper than 1/r (e.g. Li & Shu 1997; Saigo & Hanawa 1998; Tsuribe 1999; Krasnopolsky & Konigl 
2002; Shadmehri 2009). Now, we can write 



dM{r) 
dr 



27rrl]. 



The energy equation is 



7- 



1 dr 



^cl d 

— T ™w = r - A, 

r dr 



(4) 



(5) 



where F is the heating rate of the gas by dissipation processes such as turbulent viscosity and A rep- 
resents the energy loss through radiative cooling processes. The forms of the dissipation and cooling 
functions can be written as 

.dn. 



r 



(6) 



A 



1 



(7) 



7(7-1) Tcool 

where TcooI is cooling timescale. As noted in the introduction, we are interest to consider the effect of 
cooling function on the structure of self-gravitating discs. Thus, similar to Rice & Armitage (2009) we 
will study the effects of it in the case of the heating rate in disc is equal to cooling rate, F = A. 

Since fragmentation requires fast cooling, Gammie (2001) suggested the cooling timescale can be 
parameterized as /3 = VItcooI , where /3 is a free parameter Gammie (2001) showed fragmentation 
requires /3 < PctH, where pcrit ~ 3 for the adiabatic index of 7 = 2. Rice et al. (2005) performed 3D 
simulations to show the dependence of ficrit on 7: for discs with 7 = 5/3 and 7/5, PcrU « 6 — 7 and 
w 12 — 13, respectively. Recently, Cossins et al. (2010) studied /3 as a function of temperature. They 
showed that (3 has a strong dependence on the local temperature. They found that without temperature 
dependence, for radii < lOau a very large accretion rate 10^^ Af© yr~^ is required for fragmentation, 
but that this is reduced to 10^'* Mq yr~^ with cooling of dependent on temperature. So, for simplicity 
in this paper we will use a cooling timescale with a power-law dependence on temperature for study of 
the equations (l)-(5) 



Tcool - ^{^J 



■2S 



(8) 



that 5 and /?o are free parameters, and if we select Tq as a temperature of the outer part of the disc, 
then will be sound speed in there. From equation (8) and (5 = 0, we expect that ^.t^ooI becomes 
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a constant that is same with Gammie (2001) model. While non-zero 6 is qualitatively consistent with 
Cossins et al. (2010) model. It is important to stress that the above description for cooling rate is not 
meant to reproduce any specific cooling law, but is just a convenient way of exploring the role of the 
cooling timescale in the outcome of the gravitational instability. 

Here, the kinematic coefficient of viscosity can be obtained by equating of the heating and cooling 
rates 

Thus, by exploit of equation (9) we do not need to use of viscosity descriptions, such as a and /? prescrip- 
tions that are introduced by Shakura & Sunyaev (1973) and Duschel et al. (2000), respectively. Equation 
(9) implies that the kinematic coefficient of viscosity in the present model depends on physical quanti- 
ties of the system, specially cooling timescale. The kinematic coefficient of viscosity in a-prescription 
is v = acgh, where a is a free parameter and is less than unity (Shakura & Sunyaev 1973). By using 
equation (9) for a parameter we can write 

V 

CM 



1 I dil I 

\ dr\ 



1 

— -■ (10) 



7(7-1) r'^h Tcooi' 

The above equation implies that the a parameter is not a constant and varies by position and strongly 
depends on cooling timescale. We will study the a parameter in section 4 and will show that in the 
present model it increases by radii. 

As mentioned in the introduction, the gravitational stability of the disc can be investigated by 
Toomre parameter (Toomre 1964). The Toomre parameter for an epicyclic motion can be written as 

where 




2^ (12) 

is the epicyclic frequency which can be replaced by the angular frequency, VI. 

The equations of (l)-(5) and (9) provide a set of ordinary differential equations that describes phys- 
ical properties of the self-gravitating disc. Since, these equations are nonlinear, we will need suitable 
boundary conditions to solve it numerically. Thus, in next section we will try to obtain asymptotic so- 
lution in outer edge of the disc and then by exploit of this asymptotic solution as a boundary condition, 
we will able to integrate system equations inward from a point very near to the outer edge of the disc. 

Before next sections and the numerical study of the model, we shall express all quantities in units 
with values typical protostellar disc. We will choose astronomical unit {au) and the sun mass {Mq) as 
the units of length and mass, respectively. Thus, the time unit is given by ^J gu^/GMq that is equal to 
a year divided to 2tt. 

3 OUTER LIMIT 

Here, the asymptotic behavior of the system equations as r ^ i? is investigated that R is the outer radius 
of the disc. The asymptotic solutions are given by 

^w-Jtj + + ^^^^ 



Vr[r)r^-ci\j [l + a2- + ---] (14) 



6 



K. Faghei 



Of \ M^+Md^ f s \ 

^(r) - C2\j [l + az- + -- • j (15) 

CsW~C3 V ^^^R J ^ 

M{r) ^ Md,sc - / 27rr'S(r')dr' (17) 



where s = R—r, Aldisc is the disc mass, and the coefficients of c;, ai, and Sq must be determined. Using 
these solutions, from the continuity, momentum, angular momentum, energy, and viscosity equations 
[(l)-(5), and (9)], we can obtain the coefficients of q that have the following forms: 

ci = , ^ (18) 



anlil - l)/3oM(a3 - 2)(ai + aj) 
27rI]oVAf* + Mdisc{ai + aj. + 04 - 1) 

02 



C2 



1 



(19) 



C3 



where 



a37/3o(a3-2)(7-l)M 



04 = (l + a2)(l-7). 



C2 



(20) 



(21) 



The amount of mass accretion rate can be determined by observational evidences of the protoplanetary 
discs. Also, So approximately can be determined by disc mass, Mdisc ~ 7ri?^S. Thus, knowing the 
amounts of Eq and M from the observations, the value of C3 coefficient is only depended on value of C2- 
On the other hand, the value of C2 can be obtained by equation (19). Since, we have only one equation for 
coefficients of (equation 21), we will select below values for them in duration of numerical integration 
of system equations to obtain physical results 



3 3 
ai<-2+-7, 3a2=a3 = -, 04 = (1 + a2)(l - 7). 



(22) 



4 NUMERICAL RESULTS 

If the value of R is guessed, the equations by Fehlberg-Runge-Kutta fourth-fifth order method can be 
integrated inwards from a point very near to the outer edge of the disc, using the above expansions. 
Examples of such solutions for surface density, half-thickness of the disc, temperature, Toomre param- 
eter, and the viscous parameter of a as a function of radius are presented in Figs 1-5. The delineated 
quantities of T in Figs 1-4 is the mid-plane temperature and can then be determined using 




where /i = 2 is the mean molecular weight, rup is the proton mass, and kg is Boltzmann's constant. 
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Fig. 1 Surface density, thickness, temperature, and Toomre parameter of the disc as a function 
of radius, for several values of S. The surface density and the temperature are in cgs system, 
and the thickness and the distance are in au unit. The solid lines represent 6 — 0, the dashed 
lines represent 6 — 0.75, and the dotted lines represent 5 = 1.5. The input parameters are 
set to the disc mass Mdisc = O.IM©, the star mass A/, = Mq, the mass accretion rate 
M = lO~^M0yr~^, the ratio of the specific heats is set to be 7 = 5/3, and /?o = 2. 



4.1 The influences of pliysical parameters on tlie results 

The free parameters in the present model are the importance degree of temperature in cooling timescale, 
(5, the mass accretion rate, M, the parameter of /3o, the ratio of disc mass to star mass, q = Mdisc/ M^,. 

4.1.1 5 parameter 

The effects of 5 parameter on the physical quantities are presented in Fig. 1. The profiles of surface 
density and temperature show that they increase by adding 5. But, the increase of surface density is more 
than temperature. Thus, the Toomre parameter {Q cx Cs/S cx -v/T/E) decreases by adding 5 parameter. 
The profiles of Toomre parameter represent that for small 5, only outer part of the disc gravitationally 
is unstable, and the gravitational instability can extend to inner radii by adding 5 parameter. For Sent ^ 
1.5, the Toomre parameter in radii > 5 au becomes smaller than critical Toomre parameter {Qcri ^ 1) 
and the disc becomes gravitationally unstable. In the other words, the profiles of Toomre parameter 
represent the gravitational instability of the flow strongly depends on cooling timescale with temperature 
dependence. This result is qualitatively consistent with direct numerical simulations (e. g. Cossins et 
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Fig. 2 Surface density, thickness, temperature, and Toomre parameter of the disc as a function 
of radius, for several values of /Sq. The surface density and the temperature are in cgs system, 
and the thickness and the distance are in au unit. The solid lines represent /3o = 1, the dashed 
lines represent Po ~ 5.0, and the dotted lines represent (3o = 10. The input parameters are 
set to the disc mass Mdisc = O.IM0, the star mass AI^, = Mq, the mass accretion rate 
M = 10~^AfQ?/r~^, the ratio of the specific heats is set to be 7 = 5/3, and 5 — 1.0. 



al. 2010). The disc thickness increases by adding 6 parameter. It can be due to the increase of the 
temperature (ft, cx Cs cx a/T). 
Equations 8 and 9 imply that 

Since Cg > Csg the right-hand side of above equation is equal or less than one. On the other hand, 
non-zero S constrains lower viscosity for hotter regions of the disc. The study of gravitational instability 
shows that it enhances in lower viscosity (Abbassi et al. 2006; Shadmehri & Khajenabi 2006; Khajenabi 
& Shadmehri 2007). Thus, the gravitational instability can be enhanced by adding the S parameter for 
hotter regions. But there is a limitation for the amount of S parameter that we discuss it in next section. 

4.1.2 /3q parameter 

The influences of parameter of Po are shown in Fig. 2. As, we know from the simulations of self- 
gravitating disc (Gammie 2001; Rice et al. 2003), the reduce of this parameter provides conditions 
that the disc places on gravitational instability and consequently fragmentation. The profiles of surface 
density show that it does not change by adding the /?o parameter and only it shows small deviations in 
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Fig. 3 Surface density, thickness, temperature, and Toomre parameter of the disc as a func- 
tion of radius, for several values of M. The surface density and the temperature are in 
cgs system, and the thickness and the distance are in au unit. The solid lines represent 
M = 10~^Mo2/r~^, the dashed lines represent il/ = 5 x lO^^il/0yr^^, and the dotted lines 

representM = lO~^M0j/r~^. The inputparameters are set to the disc mass M^isc = O.lAf©, 
the star mass = Mq, the ratio of the specific heats is set to be 7 ~ 5/3, Pq = 10 and 
S = 1.0. 



large radii. The disc temperature increases by adding the /3o parameter. Because, the increase of this 
parameter reduces the rate of cooling. In large amount of /3o 10), the disc is gravitationally stable, 
while by reduce of its value to 5, the gravitational instability can occur in large radii, and for the small 
value of it (/3o ~ 1), we can expect gravitational instability in whole of the disc excluding near to the 
star These results are qualitatively consistent with direct numerical simulations of protoplanetary disc 
(Gammie 2001; Rice et al. 2003; Cossins et al. 2010). Also, the solutions show that the disc thickness 
increases by adding the /?o parameter. 

4.1.3 The mass accretion rate 

Rice & Armitage (2009) showed that beyond of 1 au the disc reaches a quasi-steady state in 20000 years 
and mass is redistributing itself to produce a state in which the accretion rate is largely independent of 
r. The mass accretion rate in their simulations finally reached to 10~^ — 10" A/0 /yr (see Fig 4 in their 
paper). We will study the behavior of the present model in Fig 3 for several values of the mass accretion 
rate (10"^, 5 x 10"^, and lO^^MQ/yr). The solutions imply that the disc temperature is sensitive to 
the amount of mass accretion rate and increases by adding the mass accretion rate. While, the surface 
density is not sensitive to the mass accretion rate and only shows small variations in large radii. Thus, 
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Fig. 4 Surface density, thickness, temperature, and Toomre parameter of the disc as a function 
of radius, for several values of q = Mdisc/M^,. The surface density and the temperature are 
in cgs system, and the thickness and the distance are in au unit. The solid lines represent q = 
0.05, the dashed lines represent q = 0.1, and the dotted lines represent q = 0.15. The input 
parameters are set to the star mass A/» = A/0, the mass accretion rate M ~ 10~^MQyr^^, 
the ratio of the specific heats is set to be 7 = 5/3, Po — 2 and 5 = 1.5. 



the behavior of the temperature only specifies the behavior of the Toomre parameter (Q oc -y/T/I]). 
The profiles of Toomre parameter represent that it increases by adding the mass accretion rate. Also, the 
solutions show the disc thickness increases by adding mass accretion rate, that is due to increase of the 
disc temperature. The solutions show that for a low mass accretion rate 10"'' Mq/j/f), but cooling 
timescale with temperature dependence ((5^1), the gravitational instability can occur for radii > 10 au. 

4.1.4 Mass ratio 

As noted in the introduction, semi-analytical studies of self-gravitating discs are regarding discs without 
central object. This simplification is relevant to protostellar discs at the beginning of the accretion phase, 
during which the mass of the central object is small and only self-gravity of the disk plays an important 
role. Also, this simplification can correspond to discs at large radii because the effects of the central mass 
become unimportant in the outer regions of the disc. As, the central object attends in the present model 
and its effects are not ignored. Thus, the present model does not have limitations of previous studies of 
semi-analytical self-gravitating discs and can be applied for all region of the disc. Fig. 4 represents the 
effects of the ratio of the disc mass to the star mass q = M^,/Mdisc on the present model. The solutions 
show the surface density increases and the temperature decreases. Each of the surface density increasing 
and the temperature decreasing individually can reduce the Toomre parameter Thus, we expect that 
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Fig. 5 The viscous parameter of a as a function of radius (au). The input parameters are set 
to the star mass = A/0, the mass accretion rate M = lO^^MQ-yr^^, the ratio of the 
specific heats is set to be 7 = 5/3. Left panel is for several values of Gammie's parameter 
/3o, the solid line represents /3o = 1, the dashed line represents /3o = 5, and the dotted line 
represents /3o = 10, and (5=1.5. Right panel is for several values of parameter of 5, the solid 
line represents & = 0.5, the dashed line represents 5 = 1.0, and the dotted line represents 
S= 1.5, and^o = 1-0. 



Toomre parameter decreases by adding q parameter that the profiles of Toommre parameter confirm it. 
The disc thickness profiles represent the disc thickness decreases by adding the disc mass. This property 
is qualitatively consistent with two-dimensional study of self-gravitating disc (e. g. Ghanbari & Abbassi 
2004). 

4.2 The viscous parameter of a 

In the present model, the viscous parameter of a depends on the physical quantities of the disc (Equation 
10), especially the local cooling rate which depends on the local temperature. The profiles of the viscous 
parameter of a show that it increases by radii that this property is agree with simulation results of Rice 
& Armitage (2009) and Rice et al. (2010). As, mentioned in the introduction, the minimum cooling 
timescale depends on the equation of state (Rice et al. 2005) with fragmentation occurring for TcooI < 
3ri~^ when the specific heat ratio 7 = 5/8 (Gammie 2001). Rice et al. (2005) showed that fragmentation 
occurs for a > 0.06 and this boundary is independent of the specific heat ratio 7. Left panel of Fig 5 
represents the viscous parameter of a as a function of radius for several values of the (3o parameter The 
solutions show the viscous a strongly depends on the /3o parameter. As, the a parameter decreases by 
factor of /3o- Also, the solutions for small values of /3q show the viscous a can reach to its critical value 
for fragmentation. Right panel of Fig 5 represents the viscous parameter of a as a function of radius for 
several values of the S parameter The solutions represent the a parameter excluding the outer region 
of the disc strongly depends on the S parameter In 5 ~ 0.5, the value of viscous a in whole of the 
disc is in the region for fragmentation. However, Rafikov (2005) suggested that it is extremely difficult 
to see how fragmentation can occur within 10 au even for the relatively massive discs. In S = 1.0 and 
S = 1.5, the viscous a in the inner disc (r < 10 and 40 au, respectively) is well below that required for 
fragmentation. 

The requirements for fragmentation are Q < 1 and a > 0.06 (Rice et al. 2005, 2010; Rice & 
Armitage 2009). In the present model, apparently the increase of the S parameter reduces possibility of 
fragmentation (Right panel of Fig 5). On the other hand, the increase of S parameter can place the disc 
in gravitational instability (Fig 1). Thus, by a suitable value for the S parameter, the disc can obtain two 
requirements for fragmentation. The Figs 1 and 5 imply that this value for small /?o can be between 0.5 
and 1.0. 
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5 SUMMARY AND DISCUSSION 

In this paper, we have studied self-gravitating accretion discs in presence of a Newtonian potential of 
a point mass. We have used a prescription for cooling that is introduced by Gammie (2001). But, due 
to recent results of Cossins et al. (2010), we have assumed that cooling timescale in unit of dynamical 
timescale is a power-law function of temperature. Because of, the system equations are non-linear and 
there is not self-similar solution for it. First, we have obtained asymptotic solutions for system equations 
and then by them as boundary conditions, we integrated system equations numerically. 

The solutions showed that the structure of the disc strongly depends on the present cooling function. 
As, by adding importance degree of temperature in cooling timescale, gravitational instability extends 
from outer to inner radii. The solutions showed that in the case of cooling with temperature dependence, 
the disc thickness increases. But, this change of thickness is important in region with smaller Toomre 
parameter. In the present model, the effect of physical parameters studied such as mass accretion rate, 
/3o parameter, and the ratio of the disc mass to central object mass. The results showed the structure of 
the disc is sensitive to these parameters. For example, the disc becomes gravitationally stable in larger 
mass accretion rate. While, the gravitational instability can occur in the larger disc mass. Also, the disc 
thickness increases by adding the mass accretion rate and decreases by adding the ratio of the disc mass 
to the star mass. The study of the viscous parameter a in the present model shows that it increases by 
radii that this result is consistent with direct numerical simulations (e. g. Rice & Armitage 2009; Rice et 
al. 2010). Also, the solution implies that the viscous a in the outer part of the disc becomes larger than 
its critical value (^ 0.06) that might mean condition for fragmentation. 

Here, the solutions represented that the disc thickness is very sensitive to input parameters. Thus, 
study of the present in a two dimensional approach may be interesting subject for future works. Also, it 
will be interesting to obtain a suitable & value for fragmentation by direct numerical simulations. 
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